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Abstract. We present results from a set of high (512'^ effective resolution), and ultra-high 
(1024^) SPH adiabatic cosmological simulations of cluster formation aimed at studying the in- 
■ ternal structure of the intracluster medium (ICM). We discuss the radial structure and scaling 

^SJ ' relations expected from purely gravitational collapse, and show that the choice of a particular 

halo model can have important consequences on the interpretation of observational data. The 
validity of the approximations of hydrostatic equilibrium and a polytropic equation of state are 
^ ' checked against results of our simulations. We also show the first results from an unprecedented 

^ : large-scale simulation of 500 h ^ Mpc and 2 x 512^ gas and dark matter particles. This experi- 

' ment will make possible a detailed study of the large-scale distribution of clusters as a function 

^) \ of their X-ray properties. 
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1. Introduction 

Oh| Galaxy clusters are a unique laboratory to test the hierarchical paradigm of structure 

Q i formation. They are the best probes of the large scale structure of the Universe and have 

I often been used as a diagnostic of the cosmological parameters. The intrinsic non-linear 

^ . nature of gravitational collapse and gas dynamics makes numerical simulations the most 

useful tool to study in detail the process of cluster formation and evolution. 

Simple analytical models for the structure of the ICM can be derived from the hy- 
potheses of hydrostatic equilibrium and polytropic equation of state, P (x p^- But real 
■ clusters might not be well described by these two hypothesis. For instance, kinetic energy 

I makes a significant contribution to the energy budget of merging systems, and therefore 

thermally-supported hydrostatic equilibrium ceases to be a valid approximation. Even in 
relaxed systems, this assumption is not very accurate in the outermost parts, where gas 
motions become more important. Departures from spherical symmetry can also play a 
role in the final structure of the ICM (Lee & Suto 2003) and, last but not least, there is 
no obvious physical reason for t he gas to follow a polyt ropic relation. From our numeri- 
cal experiments, we showed (see lAscasibar et al.l l|2003|) ) that hydrostatic equilibrium is 
fulfilled within ~ 20 per cent accuracy by all simulated clusters, as long as they are not 
heavily disturbed. A polytropic equation of state seems to be a good appr oximation as 
well, although its reliability near the centre is still a matter of debate (see e.g. lR,asia et al.l 
ioOS). From our data we derive a polytropic index of 7 1.18. 

We compared four different analytical halo mo dels with our si i nulat ions. The first two 
model s assume that haloes are well described by iNavarro et al.l l)l997(l and iMoore et alJ 
19991) fitt ing f ormulae, while th e other two assume that the gas follows a /3-model 



Cavaliere fc Fusco-Femianolll976|) . We consider a 'canonical' version of the /3-model, in 



which the gas is isothermal (7 = 1) and P = 2/3, and a polytropic version with 7 = 1.18 
and P = 1. The same value of the polytropic index has been used for the first two models 
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as well. Hereafter we will use the abbreviations NFW, MQGSL, BM and PBM to refer t o 
these models. For a detailed description, the reader is referred to lAscasibar et all |2£ 



2. Numerical experiments 

We have carried out a series of high-resolution gasdynamical simulations of cluster for- 
mation in a flat ACDM universe {Vl^ = 0.3; f^A = 0.7; h = 0.7; as = 0-9; = 0.02 h''^). 
Simulations have been done with the parallel GADGET code dSpringel et alJ2 00ll). with 
a nov el version of SPH in which the entropy is explicitly conserved ijSpringerfc^ernauistl 
l2002(l . In a cubic volume of 80 ft. ^ Mpc on a side, an unconstrained realization of the 
power spectrum of density fluctuations corresponding to the ACDM model was gener- 
ated for a total of 1024^^ Fourier modes. The density field was then resampled to a grid 
of 128^ particles, which were displaced from their Lagrangian positions according to the 
Zeldovich approximation up to z = 49. Their evolution until the present epoch is traced 
by means of a pure N-body simulation with 128^ dark matter particles. A sample of 
clusters selected from this preliminary low-resolution experiment were re-sim ulated with 
higher resolution by means of the multiple mass technique (see iKlvpin et a l. 2001^ for 
details). Mass resolution is then increased by using smaller masses in the Lagrangian 
volume depicted by these particles, including the additional small-scale waves from the 
ACDM power spectrum in the new initial conditions. We use 3 levels of mass refinement, 
reaching an effective resolution of 512'^ CDM particles (2.96 x 10* h^^ Mq). Gas was 
added in the highest resolved area only. The total number of particles (dark-}-SPH) in 
this area is greater than 1 — 2 x 10^ for all clusters. The gravitational softening length 
was set to £ = 2 — 5 kpc, depending on number of particles within the virial radius 
llPower et al .120031) . The minimum smoothing length for SPH was fixed to the same value 
as e. 

In order to study effects of resolution in the determination of X-ray properties of our 
clusters, we have resimulated one of the objects with 8 times more mass resolution, 
reaching an effective resolution of 1024^^ particles (i.e. rudark ~ 3 x 10^ h~^MQ\ rusph ~ 
5 X IQ^Mq). The total number of particles within the virial radius was 11, 106,465. 

The list of our simulated clusters extracted from the 80 h~^ Mpc volume span a 
relative small range in X-ray emission temperature (from 0.6 to 3 keV). In order to extend 
our numerical sample of clusters to wider temperature (mass) range, we have simulated 
a considerable much bigger volume (500 h~^ Mpc) in which a random realization of 
the ACDM power spectrum was generated with 2048^ particles. In this way, we will 
have a similar resolution for our clusters than in the previous experiments. We have 
resimulated the whole 500 h^-'^ Mpc box with different mass resolutions: 2 x 128^, 2 x 256^, 
and 2 x 512^^ dark and SPH particles. We identified all halos in the lowest resolution 
run (128'^). Using the same technique as before, we selected clusters in this run and 
resimulated them at full resolution (2048^), using 5 different species of dark particles 
and SPH with msph ~ 1-6 x 10*/i~^Mq. The clusters selected cover a range of masses 
from 2 X 10^'* — 2 x 10^^ Mq. In this regard, we could extend the temperature range of 
our cluster sample up to 11 keV. The total number of particles within virial radius in 
these clusters is comparable to the number of particles of previous simulations (~ 10^). 

We have recently finished the run with 2 x 512^ gas and dark particles that, to our 
knowledge, is one of the largest adiabatic SPH simulations of large scale structure done 
so far. The mass resolution is mdark ~ 6 x \Q^^h~^ Mq, which means that we can resolve 
from galactic halos (100-1-100 particles) to the biggest galaxy clusters (4 x 10^ particles). 
We identified a total of 4 x 10^ dark matter halos, with 10 or more particles, in this 
simulation. Well resolved halos (7000 dark particles or more) correspond to clusters with 
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Figure 1. Results from 500^ simulated volume. Left, the X-ray Temperature function and its 
observational estimates tikebe et aLii2002i'l . Middle, the M^ir —T^ relation. Right, the Lx —Tx 
relation. Stars represent the clusters resimulated with high resolution. The dotted line is a fit 
to the hottest halos ( > 6 keV). The slope is ~ 2, as in the self-similar scaling behaviour 



emission weighted temperatures > 3 keV. If we go down to 1000 particles, the halos 
resolved have temperatures > 1 keV, and temperatures go down to ^ 0.6 keV for halos 
with 500 or more particles. Due to the large simulated volume, we have a statistical 
significant sample of clusters and groups. The total number of objects with Tx > 0.6 
keV, excluding substructure, is ~ 30, 000. A total of 126 hot clusters have been found 
with Tx > 5 keV. The X-ray temperature function for clusters in this simulation is 
shown in Figure ^ together with recent observational estimates. As can be seen, the 
number density of clusters for the h ighest temperatures (> 6 keV) found in simulations 
is compatible with observations fsee lBorgani et all 1 20041) for results from non-adiabatic 
simulations). 

The resulting M — Tx and Lx — Tx relations are depicted in Figure ^ in which 
we also show 6 of the hotter clusters resimulated with high resolution. They show a 
convergence of results at least for clusters with Tx > 3 keV. The fit to the M — Tx 
relation for clusters with Tx > I keV is {M^^/Mq) = (T^j/keV)", with a = 1.56 ± 0.05 
and Mq — 4.4 x 10^^h~^AdQ. The slope is compatible with simulations which include 
non adiabatic effects, although the zero point, A/q, is a factor of 2 higher. Note however 
that we are using virial mass instead of mass for overdensity 500. This implies that the 
M-T relation is rather insensitive to the energy transfers due to non adiabatic processes 
associated to star-gas interactions. 



3. Radial structure of the ICM 

We compare the universal halo models described earlier, with results from our simu- 
lations. For NFW and MQGSL, we obtain the characteristic density and radius of each 
cluster by fitting the dark matter distribution. Gas density and temperature profiles are 
genuine predictions of these models. For BM and PBM, we fit the gas distribution and 
predict the gas temperature and the dark matter density. In Figure|21we plot the average 
dark matter and gas density profiles. Not surprisingly, both NFW and MQGSL provide 
good fits to the dark matter density. These models are able to accurately predict the 
gas density, but they are too steep at large radii due to a systematic departure from 
hydrostatic equilibrium. 

On the other hand, /3-models (BM and PBM) show a core in the dark matter density 
that is not seen in the numerical data. Moreover, they do not give an accurate description 
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Figure 2. Universal radial profiles: Upper panel dark matter and gas density. Lower panel: 
projected temperature and gas entropy. Black squares represent the numerical data, averaged 
over all clusters except major mergers. Error bars denote one-sigma scatter. Lines are used 
to plot the analytical models, and top panels quantify deviations from the simulated profile. 
Dashed lines in the projected emission- weighted temperature profile, shows the 'universal' profile 
pr oposed by L okcn ct al. (2002), while solid li ne represents the best fit to our data. Observations 
bv lDe Grandi fc Molendil II2002L circles) and IMarkevitch et all ([T99l . boxed region) are ploted 
for comparison. 



of the gas density profile. The inner regions are better described with low values of /3, 
while the outer parts require higher values for this parameter. 

One of the most striking results from our simulations is that they favour the hy- 
pothesis of a 'universal' temperature profile. As can be seen in Figure |2 the ICM is 
not isothermal, but the temperature decreases by a factor of three or four from the 
centre to the virial radius. The projected X-ray emission-weighted temperature profile 
is plotte d on the right pane l of Figure |21 We compare our simulations with previous 
work by iLoken et al.l l)2002j) based on Eulerian simulations. These authors propose a 
'universal' form Tp(r) = To(l + r/ax) * Our results are well described by this rela- 
tion. Although real clus ters seem to be consistent with a polytropic equation of state 
ijMar kevitch et al."l998') , recent observations indicate the presence of a large isothermal 
core (Dc Grandi & Molcndi 2002). Apart from this feature, which is not observed in our 
objects, adiabatic simulations arc in good agreement with observational data beyond 

0.2i?200- 

Entropy profiles are plotted on the lower panel of Figure |51 Contrary to the com- 
mon view, neither the analytical models nor the simulation data yield a pure power- 
law profile, desp ite the fact tha t purely adiabatic gasdynamics has been considered. As 
shown in Ascasibar et al.l l|2003fl . the standard implementation of the SPH algorithm can 
lead to misleading results in the inner regions due to artificial entropy losses (see e.g. 
ISuringel fc Hernauistl200^ . We find that the shape of the entropy profile does not de pend 
on th e cluster mass or temperature, in agreement with recent observations (£onman et alJ 
NFW or MQGSL models provide a better estimate of the entropy profile than the 
/3-model, but the low gas densities predicted at large radii yield a very steep slope at 
r ~ i?200- 
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Figure 3. Radial profiles for a cluster run with two different resolutions. 

4. Convergence of results 

To check for convergence of results in terms of numerical resolution, we compare, in 
Figure 01 radial profiles for one of our clusters from the ultra-high resolution simulation 
described in previous section. As it can be appreciated, the structure of the ICM is quite 
similar, at least from the virial radius down to 1% of that. The features at the most 
internal parts of the cluster are related with the different dynamical stage of the cluster. 
The positive results of this comparison is that the overall properties of the halos are well 
described with the resolution adopted in our simulations in which 500,000-1,000,000 SPH 
particles are used. 
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